Reconstruction of time-dependent coefficients: a check of approximation schemes for 
non-Markovian convolutionless dissipative generators 
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We propose a procedure to fully reconstruct the time-dependent coefficients of convolutionless non- 
Markovian dissipative generators via a finite number of experimental measurements. By combining 
a tomography based approach with a proper data sampling, our proposal allows to relate the time- 
dependent coefficients governing the dissipative evolution of a quantum system to experimentally 
accessible quantities. The proposed scheme not only provides a way to retrieve full information 
about potentially unknown dissipative coefficients but also, most valuably, can be employed as a 
reliable consistency test for the approximations involved in the theoretical derivation of a given 
non-Markovian convolutionless master equation. 
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I. INTRODUCTION 

The dissipative evolution of a quantum system inter- 
acting with an environment represents a phenomenon of 
paramount importance in quantum information science 
and beyond, as it addresses a fundamental issue in quan- 
tum theory. In general a complete microscopic descrip- 
tion of the dynamical evolution of a system coupled to 
the environment (or bath) is a complex many-body prob- 
lem which requires the solution of a potentially infinite 
number of coupled dynamical equations. According to an 
open system approach, this issue is tackled by retaining 
only basic information about the environment and de- 
scribing the system dynamics in terms of a master equa- 
tion jU [5] . The lack of a complete knowledge about the 
bath leads to master equations coefficients (MECs) which 
may be either unknown, or obtained from a microscop- 
ical derivation carried out within some approximation 
scheme. As a matter of fact it would then be highly 
appealing to devise a procedure allowing to retrieve the 
MECs starting from experimentally accessible quantities. 
This would in fact both give access to otherwise unknown 
quantities and provide a strong indication about the va- 
lidity of the adopted theoretical framework. 

So far two main dynamical regimes, Markovian and 
non-Markovian, can usually be distinguished according 
to the timescale of environment dynamics (respectively 
shorter or longer than that of the system). In [3] it 
has been shown that, in case of Markovian Gaussian 
noise, the MECs can be retrieved by means of a finite 
number of tomographic measurements by using Gaus- 
sian states as a probe. Here we want to address the 
more involved non-Markovian case. In facts, even though 



Markovian evolutions have been extensively investigated 
(see e.g. [U El SJ [5]), in general real noisy dynamics are 
far from being Markovian. Despite a growing interest 
in both theory and experiment [5J [7] a comprehensive 
theory of non-Markovian dynamics is yet to come. Ex- 
act non-Markovian master equations have been derived 
for a Brownian particle linearly coupled to a harmonic 
oscillator bath via e.g. path integral methods [H [5] or 
phase-space and Wigner function computations [10l E] . 
Analogous results have been obtained employing quan- 
tum trajectories, either exactly or in weak coupling ap- 
proximation [12H15] . In the framework of path integral 
methods, master equations have been derived both for 
initially correlated states [THIE], and for factorized ini- 
tial states in the case of weak non linear interactions |18j . 
Nevertheless, all these master equations cover only few 
cases and are not simple to solve. Indeed, it would be 
highly desirable to find an approximation scheme fully 
capturing non-Markovian features, as in general different 
approximations may lead to irreconcilable dynamics pQ. 

When dealing with a time-dependent generator one can 
face two, in principle distinct, kinds of problems. On one 
hand, taking for granted the functional form of the dis- 
sipative generator, one might be interested in retrieving 
the time-independent parameters which characterize it. 
This issue has been tackled in [T^ , were it has been pre- 
sented a first generalization of the tomography-based ap- 
proach proposed in [3] to time-dependent dissipative gen- 
erators. There it has been shown that, under the assump- 
tion of Gaussian noise, the time-independent parameters 
(TIPs) characterizing the time-dependent MECs (whose 
functional form is previously known) as for example the 
system-bath coupling, temperature and bath frequency 



cut-off, can be obtained with a finite set of measurements. 
On the other hand, one might be interested in the more 
general problem of reconstructing the functional form of 
the dissipative generator itself. This might be the case ei- 
ther because the time-dependence of the dissipative gen- 
erator is completely unknown or because one wants to 
test the validity of the theoretical assumptions at the ba- 
sis of the dissipative model. 

In this paper we propose an experimentally feasible 
procedure which allows the full reconstruction of the 
time-dependence of convolutionless non-Markovian gen- 
erators. Indeed, whenever the assumption of Gaussian 
noise is satisfied, by combining a tomography based ap- 
proach and a proper choice of data sampling, start- 
ing from a finite and discrete set of measurements it 
is possible to obtain global information about the time- 
dependent MECs. We envisage two main applications of 
the proposed scheme: both a full reconstruction starting 
from previously unknown MECs, and a sound, reliable 
and complete consistency test of the theoretical assump- 
tions made in deriving the dissipative generator. 

The paper is organized as follows. In section [TT] we 
review the tomography based approach introduced in [3] 
allowing to reconstruct the first two cumulants of a Gaus- 
sian state at any fixed time with a finite amount of mea- 
surements. In section Unl we introduce the main lines of 
our reconstruction scheme for two in principle distinct 
cases based on the purpose of the reconstruction and on 
the available prior knowledge of the MECs. In section 
|IV| two different implementations, integral and differen- 
tial, of the reconstruction procedure are presented and 
their application is discussed. To provide an explicit ex- 
ample of both approaches, in section [V] they are applied 
to the specific model of a brownian particle interacting 
with an Ohmic bath of harmonic oscillators. In section 
|VI| we summarize and discuss our results, and finally in 
Appendix[A]we provide some more details about the sam- 
pling theorems involved in our reconstruction scheme. 
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II. THE T-C PROCEDURE 



In this section we briefly review the main steps of the 
procedure introduced in [3] which allows to relate the 
evolved first and second order cumulants of a Gaussian 
state to tomographic measurements. We consider a mas- 
ter equation with unknown coefficients which generates 
a Gaussian Shape Preserving (hereafter GSP) dissipative 
evolution or, in other words, a Gaussian map. We then 
focus on the time evolution of the first and second cumu- 
lants of a Gaussian state, which completely determine 
its dynamics. Given a master equation governing the 
evolution of the density matrix of the system p(t), the 
dynamical equations for the first two cumulants formally 
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The explicit form of these equations depends on the 
adopted master equation whose coefficients, in general, 
are either unknown or derived by means of phenomeno- 
logical assumptions. Inserting into Eqs. (1) the explicit 
expression of p(t) and then solving for the MECs, allows 
to write these coefficients (or the differential equations 
they satisfy) as a function of the first two cumulants 
of a Gaussian probe. Hence, by measuring these two 
quantities, one can gain indirect experimental informa- 
tion about the MECs. To this aim, symplectic quantum 
tomography arises naturally as a key tool to perform the 
required measurements. In particular we will use the 
symplectic transform, or M 2 -transform ;20J, that finds 
its natural implementation in experiments with massive 
particles, as for example those to detect the longitudi- 
nal motion of neutrons [21] and to reconstruct the trans- 
verse motional states of helium atoms [22] • We also note 
that the symplectic transform is equivalent to the Radon 
transform |23) which is experimentally implemented by 
homodyne detection in the context of quantum optics 
|24) . We now recall a procedure introduced in [3] allow- 
ing to connect, at any fixed time t, tomographic measure- 
ments and the first two cumulants of a Gaussian probe 
by means of a small number of detections. 

We begin by switching from the space of quantum 
states, i.e. the Hilbert space, to phase-space, which can 
be done by means of the Wigner map. In facts, a quan- 
tum state p(t) can be represented on phase-space in terms 
of its Wigner function W defined as 
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The Wigner function of a Gaussian state is itself a Gaus- 
sian function on phase space. Hence, if the dissipative 
evolution is GSP, at any time the Wigner function of the 
evolved Gaussian state can be written as 
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In general performing a tomographic map then amounts 
to project the above Wigner function along a line in 



phase-space 

X- m -vp = 0, (4) 

or, in other words, to compute its symplectic transform 

w(X, p, v) = (S (X — pq — up)) 

W(q,p, t)S (X — pq — i/p) dqdp. 



(5) 

For the case of Gaussian states, the symplectic transform 
is again a Gaussian function: 



w (X, fi, v) = 
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where the second cumulants always obey the constrain 
Aq 2 [i 2 + Ap 2 v 2 + 2a(q,p) t pv > as a consequence of the 
Schrodinger- Robertson relation [25] . which represents a 
generalization of the Heisenberg principle. At this point 
one could wonder whether, due to experimental errors, 
a violation of the uncertainty principle might be ob- 
served. This may happen if measurements are performed 
on states almost saturating the inequality, i.e. on pure 
states. In our case measurements are performed on states 
undergoing a dissipative evolution which typically are far 
from being pure hence from saturating the uncertainty 
relation. Furthermore any additional noise of statistical 
origin will have the effect of moving the reconstructed 
state further away from the boundary, as noted in [26] . 

The symplectic transform (pSJ) provides a relation be- 
tween the first and second order cumulants and the tomo- 
gram values along each line. In facts, performing a tomo- 
graphic measurement consists in choosing a pair (fl,v), 
i.e. a line in phase space, and in computing m{X, p 1 v) 
for a given value of X. In J3] it has been shown that 
by choosing the lines corresponding to position and mo- 
mentum probability distributions (i.e. {fi,v) = (1,0) 
and (0,1)) then, for any fixed t, at most four points 
along each line (i.e. tomogram) are needed to retrieve 
the first and second cumulant of the associated variable. 
Analogously, the covariance of the two variables is ob- 
tained by measuring at most two points along the line 
(/2, u) — (l/v2) l/v2)- In overall, starting from a GSP 
master equation, the first and second order evolved mo- 
menta at time t can be obtained via a total amount of 
eight or at most ten points measured in three different di- 
rections. In the following, we will refer to this procedure 
as the tomograms-cumulants (T-C) procedure. 

To summarize, at each time t, the T-C procedure com- 
bined with Eqs. ([I]) provides the sought-for bridge be- 
tween dynamical parameters and measurable quantities. 
We note that, in the Markovian case, this combination 
is straightforwardly performed by inverting Eqs. (IT]) 3J. 



Finally, we emphasize that even though the assumption 
of Gaussian noise might be seen as an idealization, it is 
actually well fitted for a significant number of models 
[TJ [2] and small deviations from Gaussianity would only 
introduce small and controllable errors. Also, Gaussian 
probes are quite straightforward to produce either with a 
laser or an ordinary light source (obtaining, respectively, 
a poissonian or a thermal distribution) [271 . 



III. SKETCH OF THE RECONSTRUCTION 

The T-C procedure allows to gain information about 
the unknown MECs -via experimentally accessible 
quantities- locally in time. This is enough if the GSP 
dissipative generator is time-independent, as the infor- 
mation needed to fully reconstruct the unknown MECs 
is obtainable by performing tomographic measurements 
at one arbitrary instant of time, as shown in [3]. This is 
no longer true when facing a time-dependent generator. 
In this case to fully reconstruct the unknown MECs we 
in principle need to gather information globally in time. 
Now, taking the T-C procedure as a starting point, the 
problem we tackle is how to retrieve the full information 
we need by combining measurements performed on a fi- 
nite and discrete set of times, i.e. starting from partial 
information. 



A. the dissipative generator 

We focus on the reconstruction of the time-dependent 
MECs of the following class of master equations 
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where the unknown MECs are A(t), D qq {t), D pp (t), 
D qp (t). The system Hamiltonian H is chosen to be at 
most a second-order polynomial in the position and mo- 
mentum operators p, q |28j 



H = H + -(qp + pq) 



1 



2m 



ft2 



(8) 



where the time dependence which may be introduced by 
the Lamb shift term has been neglected. This is typi- 
cally justified as most of the times either the Lamb shift 
is negligible or the Hamiltonian part of the dissipative 
generator reaches its asymptotic value on a much shorter 
timescale compared to the non unitary part [T]. In over- 
all, the choice of operators in Eq. (l7| represents a natural 



generalization of the GSP time-independent master equa- 
tion introduced in [28 . The investigation of this time- 
dependent class of master equations is further motivated 
by the existence of a wide range of models obeying a GSP 
dissipative dynamics of this form [9"HTTl [T5HT7] . 

As a side remark we note that there may be some am- 
biguity in literature about the relation between time- 
dependent dissipative generators and Markovian/non- 
Markovian dynamics. In facts, some authors classify as 
non-Markovian only generators containing a convolution 
integral. It has recently been proved in 29 that these 
generators can be mapped into convolutionless ones. The 
distinctive feature of non-Markovianity becomes then the 
dependence of the convolutionless generator on t — to 
where io is the initial time. According to this approach 
a time-dependent generator of the kind ([7]) could be con- 
sidered as Markovian. However, following a consistent 
part of literature, e.g. [HI EH HH UHl HH HH ED] , we will 
term non-Markovian also convolutionless time-dependent 
generators as the one in Eq. §f§. 



B. the reconstruction scheme 

Our strategy towards the reconstruction of the MECs 
in Eq. ufh is made up of three main steps: 

1. use the T-C procedure to get indirect measure- 
ments of the evolved cumulants at different times; 

2. use these measurements to retrieve the values of 
the MECs (or functions of them) at those times 
exploiting the connection between MECs and first 
cumulants of a Gaussian probe, obtained by em- 
ploying the dynamical equations (flj); 

3. starting from the obtained discrete and finite set of 
values, reconstruct the full expression for the MECs 
by applying proper sampling theorems. 

In particular, useful for our purposes will be the Nyquist- 
Shannon theorem (3 11 I32j and one of its more sophisti- 
cated generalizations involving an additive random sam- 
pling [311 1331 I34j . In principle we can distinguish two 
different applications of the reconstruction procedure: 
check of the a priori assumed time-dependence of MECs 
(Case I) or complete reconstruction of MECs with no a 
priori assumptions (Case II). 

Case I: we assume a priori a certain time dependence 
of the MECs as a consequence, for example, of a micro- 
scopical derivation of the master equation. In this per- 
spective, we are interested in experimentally reconstruct- 
ing the MECs to check the validity of the approximations 
made. A mismatch between the assumed and the mea- 
sured MECs would in fact provide a strong evidence of 
the breakdown of the adopted approximation scheme. In 
this case a full knowledge of the MECs is assumed, includ- 
ing that of the TIPs involved, such as the bath frequency 
cut-off, system-bath coupling, etc. The TIPs can be ei- 
ther assumed or previously reconstructed [19] . In this 



case, given the prior knowledge of the bandwidth associ- 
ated to the function to be reconstructed (i.e. the width 
of its Fourier transform), the suitable sampling theorem 
can be chosen accordingly. If the function is band-limited 
then to obtain an exact reconstruction it is enough to 
apply the simplest sampling theorem, i.e. the Nyquist- 
Shannon theorem (see Appendix [A| . The function can 
hence be reconstructed starting from a discrete set of 
values spaced according to the width of its Fourier trans- 
form. If the bandwidth is not limited one could truncate 
it and still apply the same procedure, which would then 
be affected by the so-called aliasing error. To minimize 
it one can perform a proper truncation. Alternatively a 
more general additive random sampling, that avoids the 
aliasing error (see Appendix [A| , can be employed. 

Case II: here we want to fully reconstruct the MECs, 
or derived functions, with no previous assumption on the 
dynamics, i.e. the MECs are fully unknown. In general, 
this implies no prior knowledge of the bandwidth asso- 
ciated to the function to be reconstructed. In this case 
we must resort to an additive random sampling (see Ap- 
pendix [A|. If on one hand this procedure involves func- 
tion averages with respect to the probability of draw- 
ing n sampling times (i.e. more involved measurements), 
on the other it does not require prior knowledge of the 
bandwidth and is an alias-free sampling. To obtain the 
averages of the function we should in principle perform 
measurements over a continuous interval of time, as the 
reconstruction is proposed with continuous random pro- 
cesses. In practice, every experimental apparatus em- 
ployed to record and process the data has a dead work- 
ing time, such that the random process will be discrete in 
time, no matter how dense, thus introducing an intrinsic 
source of error in the procedure. 

As a conclusive remark, we note that in both cases the 
set of measurements required turns out to be discrete 
but in principle infinite, as the reconstruction should be 
performed over the whole real axis. This number can be 
made finite by invoking the largely reasonable physical 
condition of a finite observation time. 



IV. DIFFERENT PROCEDURES 

In this section we illustrate in more detail how to im- 
plement the above sketched reconstruction scheme. The 
starting point are the general equations (fTl) and (I7j which 
can be expressed in compact matrix form as 



dt 
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S(t) = (M-X(t)I 2 )S(t), 



(9) 



X(t) = (R-2X(t)I 3 )X(t) + D(t), (10) 



where / 2 (3) is the 2(3)-dimensional identity matrix. The 
vectors S(t) and X(t) correspond, respectively, to the 



first and second order cumulants 
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and the matrices M and R contain the Hamiltonian pa- 
rameters 
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Finally, D(t) is the diffusion vector 

(mu)D qq (t) s 
mui 
D qP (t) 



(13) 



such that the MECs to be reconstructed are \(t) and 
D(t). The dynamical evolution of the first cumulants, 
Eq. (p)]), only depends on the friction coefficient A(i) 
whereas that of the second cumulants, Eq. (10), depends 



on the whole set of MECs. At this point we can distin- 
guish two different approaches towards the reconstruc- 
tion, i.e. integral and differential. 



A. Integral approach 



where A(t,t') = J t ,dt"X(t"). Again, considering A(t) 
as a known quantity from Eq. (|14[), in Eq. (18) exper- 



imentally accessible quantities and MECs are grouped 
on different sides, respectively right and left-hand. The 
right-hand sides of both equations can thus be regarded 
as experimental measurements of the corresponding left- 
hand sides. The first step is then to compute the left 



hand side of Eq. (14) using the friction coefficient X(t) 
provided by the assumed model. A theoretical value for 
A(i), say A theor (t), is obtained and its Fourier transform 
is performed. Once obtained the bandwidth associated to 
A* r (£) and consequently the amount of points in time 
required for the full reconstruction (according to the cho- 
sen sampling theorem) the following step is to evaluate 
S(i) (Eq. (15)) on these points by successive applications 
of the T-C procedure. Finally, the last step is to apply 
the sampling theorem to the experimental data and to re- 
construct the left-hand side of Eq. (fl4J) , let us denote it 
by A ex P\t). If the match between A Wor (i) and A expt (t) 
is positive within the desired accuracy, one can then pro- 
ceed further to the reconstruction of the left-hand side 
of Eq. (18) according to the same procedure. We note 



that the integral approach is also feasible in case of time- 
dependent Hamiltonian parameters (i.e. m(t), ui(t), S(t)) 
as long as the generator remains GSP. In general this 
kind of generators requires a numerical evaluation of the 
integrals in Eqs. d9| and (10), according to which the 
suitable sampling theorem must be chosen. 



The formal solution of Eq. p)]) can be written as 

'^(0)\ 
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where the suffix j — 1,2 labels the two components of 
the vector 



S(t) 
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We note that in Eq. ([14]) the measurable quantities and 
the unknown MEC appear on different sides. Analo- 
gously, we rewrite Eq. ( 10 1 as 



±X(t) = D(t), 



where 
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Eqs. (15) and (17) are always invertible, provided one 
sets the quantity \/S 2 — uS 2 = r\ to the val ue i Q whenever 
■q 2 < [28j . The formal solution of Eq. (16) is given by 



X(t) = X(Q) + /„ dt'D(t'), which, using Eq. (J17j), can be 
recast in terms of X(t) as 

"* dt'e- 2A ^^he^-^ R D(t') 

= X{t) 
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B. Differential approach 



The T-C procedure allows to measure not only the 
cumulants of a given Gaussian state but also their first 
time derivatives. Indeed, by measuring each cumulant 
at two different times t and t + St, its derivative can be 
estimated via its finite incremental ratio. For instance 
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where the amount of time St is taken as the smallest time 
interval which can elapse between two different measure- 
ments. Once the cumulants and their derivatives at fixed 
times are given as experimental inputs, Eqs.Q and (10) 
become linear algebraic equations whose unknown are the 
MECs. Hence from Eqs. ^ and ( 10 ) we obtain the de- 



sired experimentally accessible estimations for the MECs 
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The number of points required to fully reconstruct the 
MECs depends on the sampling theorem employed. In 
Case I, given the prior knowledge of the bandwidth, if 
the function is band limited or it can be truncated on an 
effective compact support outside of which the contribu- 
tions are negligible, the Nyquist-Shannon theorem can 
be applied. Otherwise one must resort to the additive 
random sampling theorem, which must be in general 
employed in Case II. 

As the integral approach, the differential approach is 
also suitable in case of time-dependent Hamiltonian 
parameters as long as the generator remains GSP. 
The Hamiltonian parameters would then be included 
among the MEPs to be reconstructed, thus raising 
their number to seven. Since Eqs. (pi and (10) are five, 
two more equations would be required as for example 
those for two higher order cumulants, e.g. Tr(p(t)q 3 ) 
and Tr(p(i)p 3 ). This would not increase the number 
of experimental measurements, since the higher order 
moments and cumulants of Gaussian states are com- 
pletely determined by the first and the second cumulants. 



C. Comparison between the two approaches 

Let us now briefly compare the two procedures de- 
scribed in this section. The differential approach requires 
more experimental measurements compared to the inte- 
gral approach, whereas the latter might result more in- 
volved from a computational point of view. For example, 
within the frame of Case I, the computation of the first 
member of Eq. ( 18 1 might require some numerical or ana- 
lytical approximations. As any kind of approximation in 



principle reduces the accuracy of the reconstruction, the 
differential approach would be a better choice. On the 
other hand, if the computation of the integral functions 



the integral procedure should be preferred, since it re- 
quires a lower number of interactions with the physical 
system. 

In Case II, the integral approach could be employed to 
reconstruct the left hand sides of Eqs. (14)-(18), which 



are functionals of the unknown MECs. By means of time 
derivatives and linear operations on the reconstructed 
functions, the MECs can finally be retrieved. However, 
the time derivatives may amplify the error of the recon- 
struction. For instance, the Nyquist-Shannon theorem 
requires a truncation of the Fourier's frequencies, thus in- 
ducing an oscillating behavior of the reconstructed func- 
tions, i. e. introducing the so-called aliasing error (see 
Appendix A). Even if the oscillations around the mean 
(true) functions are small, the time derivative may in- 
crease them. Hence, either one performs a better recon- 
struction (e.g. a larger truncation or a random additive 
sampling requiring a larger number of measurements) or 
one adopts the differential approach, thus directly recon- 
structing the MECs. 

In general, one could say that the integral procedure is 
more advantageous in terms of number of measurements, 
but requires the ability of solving potentially involved an- 
alytical expressions. The differential approach, instead, 
is more advantageous from the point of view of versatility, 
as it allows to deal in a straightforward way with com- 
plex generators (i.e. exhibiting time-dependent Hamil- 
tonian parameters) , at the expenses of a higher number 
of measurements. Therefore, which of the two proposed 
approaches proves better, strictly depends on the specific 
case under investigation. 



V. AN EXAMPLE 

To provide an example of how to implement the pro- 
posed integral and differential procedures, in this sec- 
tion we apply both to a specific model i.e. within the 
perspective of Case I. As a model to be tested we con- 
sider a harmonic oscillator of frequency ui (our system- 
particle) linearly coupled to an Ohmic bath of harmonic 
oscillators. The TIPs of this model are the system-bath 
coupling constant a, the Lorentz-Drude cut-off lo c and 
the temperature T [1 j. Starting from a superoperatorial 
version of the Hu-Paz-Zhang master equation [5], in the 
weak-coupling limit (i.e. up to a 2 ) and using the rotat- 
ing wave approximation (i.e. an average over the rapidly 
oscillating terms) , for this model a generator of the form 
of Eq. ([7]) can be obtained [JUl ES]- More specifically, 
the Hamiltonian parameter S is set to and the MECs 
provided by the model are 
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in Eqs. ( 14 )-( 18 1 does not exhibit remarkable difficulties, 
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where at high temperature T 
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With this choice of MECs the dissipative generator is 
GSP [30]. The Markovian limit is recovered when the 
time dependent parameters X theor (t) and A theor (t) reach 
their stationary values, i.e. at times larger than l/w c 
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Usually, when studying quantum Brownian motion, one 
assumes ui c /ui >■ 1 with ui c — > oo, corresponding to a 
natural Markovian reservoir. In this limit, the thermal- 
ization time [30j is inversely proportional to the coupling 
strength, while for an out-of-resonance engineered reser- 
voir with uj c /lj <C 1 (i.e. highly non Markovian), the 
thermalization process is slowed down. 

Microscopic derivations of Master Equations usually 
give a time dependent renormalization of Hamiltonian 
parameters (to, uj, S). However, they are negligible for 
this benchmark model in the considered limits. 



A. Integral approach 

Fo llow ing the steps of the inte gral approach we solve 
Eq. (HI) using \ theor (t) (Eq. @) and we get 
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The obtained function does not belong to the functional 
space L 2 ((R)) thus not matching the condition for the 
applicability of the Nyquist-Shannon theorem (see Ap- 
pendix AJ). However, as we are interested in a finite time 
interval, we can restrict the support of A theor (t) to [0, i] 
and define: 
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The function A theor (t) is in L2OR) thus can be recon- 
structed by the Nyquist-Shannon theorem. The discon- 
tinuity at t, by inducing the Gibbs phenomenon (i. c. 
a finite Fourier sum overshoots at the jump), might at 
this point constitute a source of error. This problem can 
be anyway overcome by slightly restricting the domain in 
which the reconstruction of A theor (t) can be trusted to 
[0, i — £] with £ > 0. Performing the Fourier transform 
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FIG. 1: (Color online) Comparison A expt (t) (dotted) and i 
A or (t) (continuous), in the Markovian regime, i.e. ui c /lo = 
10. The support of A theor {t) is [0, 12/w c ] and the reconstruc- 
tion is trusted in [0, 10/w c ], with a — 0.1 and T = 10{fuj)/kB- 
In Fig. 1(a) contributions smaller than 0.1% of the maxi- 



mum value of the Fourier transform are neglected (bandwidth 
W — 19.4/2-7T, 7 reconstruction points) while in Fig. 1(b) 
those below 0.01% (bandwidth W = 196/27T, 74 reconstruc 
tion points). 



of A theor (t) we obtain 
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From the above equation we note that A theor (t) is not 
band limited, i.e. F[A theor ]{s) is not defined on a com- 
pact support. Nevertheless, since it is symmetrically de- 
creasing around s — (modulo a negligible oscillating 
behavior for higher frequencies), we can define an effec- 
tive limited bandwidth and apply the Nyquist-Shannon 
theorem by truncating its support to a symmetric inter- 
val around zero. In Figs. Ill and pi we compare A theor (t) 
and A ex pt (t) (obtained by applying the reconstruction 
formula (Al)) in both the Markovian (u c /ui — 10) and 
the non-Markovian (uj c /uj = 0.1) limit. To respect the 
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FIG. 2: (Color online) Comparison between A expt (t) (dot- 
ted) and A theor (t) (continuous) in the non-Markovian regime, 
0.1. The support of A theor (£) is [0, 12/w c ], with 
and the reconstruction is trusted in [0, 10/oj c ], 

contributions 
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with a = 0.1 and T = lQ(tko)/k B 
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In Fig 

smaller than 0.1% of the maximum value of the Fourier trans- 
form are neglected (bandwidth W = 0.16/2-7T, 6 reconstruc- 
tion points) while in Fig. 2(b) those below 0.01% (bandwidth 

w 



1.66/27T, 64 reconstruction points) 



conditions of weak coupling and high temperature, in 
both plots a is set to 0.1 and T to 10(Hui)/kB. The re- 
construction of A expt (t) is trusted in [0, i— £ = 10/w c ]. As 
already mentioned, the truncation of the Fourier trans- 
form required by the Nyquist-Shannon theorem is the 
source of the so-called aliasing error in the reconstruc- 
tion. In order to provide evidence of this phenomenon 
we compare the results obtained for two different effec- 
tive supports of the Fourier transform 



in Figs. 1(a) 
and |2(a)| we consider only contributions higher than 0.1% 
of the maximum value, whereas in Figs. l(b)| and 2(b) 
those higher than the 0.01%. As expected, the second 
case returns a better approximation and the aliasing er- 
ror becomes negligible, thus not requiring the application 
of additive random sampling. For more involved Fourier 
transforms, a more general truncation criterion would be 
taking as effective support the region where the integral 
of the Fourier transform is greater than a chosen thresh- 
old value. We skip the reconstruction of the remaining 
MECs as the procedure is analogous, the only difference 
being that the integrals on the left-hand side of Eq. ( 18 ) 
require a numerical evaluation. 



B. Differential procedure 

Let us now apply the differential procedure to our 
model. We note that in this specific example, as the four 
MECs to be reconstructed effectively reduce to two (A(£), 
A(£)), we could in principle allow for time dependent 
Hamiltonian parameters (rn(t),u(t)) without increasing 
the number of equations. However, to be consistent with 
the previous section and with Ref. |30j . we will con- 
sider time-independent Hamiltonian parameters. For our 
model Eqs.(20) and (23 1 read 
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(31) 



Once again we restrict ourselves to a compact support 
[0, i\ (trusting the reconstruction in [0,i— £]) by defining 
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The Fourier transforms of the restricted functions read: 
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As in the previous section, the supports of both 
T[X expt ](s) and T[A expt ](s) are not compact, such that 
we need to truncate them. Having already recon- 
structed A expt (t) in the previous section, we focus here 
on the remaining MEC. In Figs. [3] (Markovian regime) 
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FIG. 3: (Color online) Comparison between A expt (t) (dot- 
ted) and A theor (t) (continuous) in the Markovian regime, i.e. 
uj c /uj = 10. The support of A expt (t) is [0, 12/w c ] and the 
trusted interval is [0, 10/w c ] with a — 0.1 and T = 10(fko)/kB- 
In Fig. 3(a)| contributions smaller than 0.1% of the maxi- 
mum value of the Fourier transform are neglected (bandwidth 
W = 19.5/2-7T, 7 reconstruction points) while in Fig. 3(b) 
those below 0.01% (bandwidth W = 73/2-7T, 34 reconstruc 
tion points). 
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FIG. 4: (Color online) Comparison between the recon- 
structed function A expt (t) (dotted) and its theoretical be- 
havior A i '' e ° r (i) (continuous) in non-Markovian regime, i.e. 
ui c /uj = 0.1. The support of A expt (t) is [0, 12/w c ] and the 
trusted interval is [0, 10/w c ] with a = 0.1 and T — lO(huj)/kB- 
contributions smaller than 0.1% of the maxi- 



In Fig. 4(a 



mum value of the Fourier transform are neglected (bandwidth 
W — 1.32/27T, 50 reconstruction points) while in Fig. 4(b) 



those below 0.01% (bandwidth W = 3/2tt, 114 reconstruction 
points). 



and H (non-Markovian regime) the reconstructed func- 
tion, A theor (t), (obtained by applying the reconstruc- 
tion formula Eq. ( |Al[ )) is compared with the theoret- 
ical curve A theor (t) (Eq. (25)) within the time interval 
[0,i — £ = 10/Wc]. As in the previous section, the cou- 
pling constant and the temperature are set to a = 0.1 and 
T = 10(fku)/kB- Again, in order to provide an indication 
of the sensitivity of the differential procedure to different 
truncations of the Fourier transform (aliasing error), in 
Figs. |3(a)| and |4(a)| the contributions smaller than 0.1% 
of its maximum value have been neglected, whereas in 
Figs. |3(b)| and [4(b)1 those smaller than the 0.01%. Anal- 
ogously to what happened for the integral procedure, the 
truncation to 0.01% returns a better reconstruction and 
the aliasing error becomes almost negligible (thus not re- 
quiring additive random sampling). 



VI. CONCLUSIONS 

In this paper we have proposed a procedure to re- 
construct the master equation coefficients in a non- 
Markovian scenario by means of a finite and discrete 



set of tomographic measurements. In particular, our 
proposal applies to non-Markovian convolutionless Gaus- 
sian Shape Preserving dissipative generators. While the 
tomogram-to-cumulants procedure provides the local-in- 
timc relation between the evolved cumulants of a Gaus- 
sian probe and experimental measures, a careful inter- 
play between dynamical equations and sampling theo- 
rems tells us how these measurements must be repeated 
in time to access full global information about the time- 
dependent master equation coefficients. We have then 
proposed two different approaches, integral and differen- 
tial, towards full reconstruction, and provided an explicit 
example of their application. According to whether or 
not one has previous knowledge of the master equation 
coefficients to be reconstructed, our proposal can be cast 
within two in principle distinct perspectives. If on one 
hand it provides a tool to gain information about fully un- 
known quantities, on the other, if any a priori knowledge 
is available, it represents a consistency test of the theo- 
retical assumptions associated to a given model, there- 
fore providing reliable information about the overlap be- 
tween the real and the assumed dissipative dynamics. In 
case one has previous knowledge of the MECs and can 
perform the intergrals in Eqs. (14 1,(18 1, the integral ap- 
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proach proves more efScient in terms of the amount of re- 
quired measurements, otherwise the differential approach 
represents the most appropriate one even though it re- 
quires more measurements. Furthermore, along this line 
of thought our procedure might also prove useful to test 
the theoretical predictions associated to a given model, 
such as the crossing of the Lindblad-non Lindblad bor- 
der investigated in |30j . Our proposal opens up several 
interesting questions which are going to be the subject of 
further future investigation. In facts how our approach 
can be recast within an estimation theory perspective, 
along the lines of [33], represents a relevant open scenario. 
Another relevant point to investigate is whether the pro- 
posed protocol can be enhanced by employing entangled 
Gaussian states as a probe. Finally, whether or not our 
proposal can be further generalized and employed in pres- 
ence of memory kernels represents a non trivial question. 
Indeed, finding a way to reconstruct with a finite amount 
of measurements the coefficients of non-Markovian dissi- 
pative generators with memory represents both a highly 
challenging and interesting task which deserves further 
future investigation. 
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Appendix A: Sampling theorems 

Here we provide some details about the Nyquist- 
Shannon and the additive random sampling theorems, 
for a full review on this topic see [3_TJ [32] . The simpler 
sampling theorem, known as Nyquist-Shannon theorem, 
deals with functions whose Fourier transform has a com- 
pact support. The theorem allows these functions to be 
reconstructed starting from a discrete and infinite set of 
values. The theorem is stated as follows: 
Theorem 1. (Nyquist-Shannon theorem) If a function 
F(t) has no frequency higher than W , that is the support 
of its Fourier transform is contained in [— 2ttW, 2ttW], 
then it is completely determined by giving its ordinates 
at a series of points spaced -^m apart. 

The reconstruction formula is 
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The most frequent sources of error are the truncation 
error, the aliasing error, the round-off error and the 
jittering error. The truncation error arises from con- 
sidering a finite sampling instead of an infinite one, 
as in practice it is unfeasible to sample and store an 



infinite number of values (unless some regular behavior 
of the function to reconstruct can be postulated). The 
aliasing error occurs whenever a function which is not 
band-limited is reconstructed by means of procedures 
suitable for band-limited functions which is often the 
case since band-limited functions are very peculiar. 
Whenever a truncation or an aliasing error occurs, it 
means that there are different functions matching the 
exploited sampling, anyway there exist conditions to 
bind this kind of errors [3TJ [32] ■ Whereas the two pre- 
vious errors are sampling-related, whereas the round-off 
and the jittering error are linked to the precision of the 
experimental apparatus. In particular, the round-off 
error is caused by errors affecting the sampling values 
and, finally, the jittering error is due to errors affecting 
the sampling times. 

A remarkable generalization of the sampling theorem 
is developed in [33[ I34j and requires an additive ran- 
dom sampling. The function is sampled at additively 
randomly chosen points t n = £ n _i + 7«,, where {jn} is 
a family of independent identically distributed random 
variables, whose probability distribution p(t) obeys the 
following constraints 

p(t)eL 2 (R), p(t)=0 for t < 0, E[ 7 „] = h< oo, 

(A2) 
where E[-] is the expected value and h is the average spac- 
ing of the sampling. Additive random samplings may 
be alias-free, namely they may provide reconstructions 
without any aliasing error. As a consequence, non band- 
limited functions may as well be reconstructed. Exam- 
ples of alias- free samplings are discussed in [33l [34] . The 
following theorem provides the condition for an additive 
sampling to be alias-free: 

Theorem 2. An additive random sampling is alias- free 
if the characteristic function 4>{lo) = E[e lut ] takes no val- 
ues more than once on the real axis. 
Conversely, if the characteristic function 4>(uj) takes the 
same value at two different points of the open upper half- 
plane, then aliasing occurs with an additive random sam- 
pling. 

The reconstruction formula is a bit involved and is 
given in [33j [34] . We emphasize that in this case the in- 
put sampling is not a countable set of values but rather 
a countable set of averages as the sample spacing is a 
random variable itself and we make an average over all 
possible spacings. If the spacing distribution is highly 
picked on h, i.e. p(t) = 5(t~h), the spacing variance van- 
ishes E[7 T 2 J = E 2 [7„], and the reconstruction procedure 
recovers the usual one with equally spaced samplings. If 
the spacing variance is comparable with its mean value 
E[7 r 2 J — E 2 [7 n ] ~ h? , then all spacings are highly prob- 
able and we get an effective continuous sampling. An 
unavoidable source of error in the implementation of this 
sampling theorem is that it involves functions defined on 
arbitrarily small times [331 134] whereas any experimen- 
tal apparatus exhibits a dead working time interval to 
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record and process data. Finally we note that we want 
to reconstruct functions defined only for t > hence we 
do not encounter any lower truncation error. The upper 
truncation error can be also avoided by reconstructing 
slightly different functions: F(t)(l — 9(t — £)) instead of 
F(t). The differences arise from the integral transforms 
involved in the reconstruction formula, but if we are in- 
terested in reconstructing functions in the experimentally 



accessible timescales, values greater than the threshold t 
of experimentally reachable times can be neglected. The 
same trick can be exploited to avoid errors due to a fi- 
nite input size for an additive random sampling. Indeed 
if the sampling probability is well localized, each p n (t) 
is localized as well. The probabilities p n (t) localized at 
times larger than the threshold t do not contribute and 
the corresponding /„ vanish. 
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